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1. Introduction 

Yang-Mills theories are the cornerstone of the standard model. These theories are formulated 
in terms of a highly redundant set of variables, due to the local gauge symmetry of the action. In 
case of Yang-Mills theories with matter, the variables are the gauge fields (x) and the matter 
fields Q{x). The calculation of physical observables is hampered by the fact that these fields are 
not directly related to physical particles. This can be easily seen from the fact that they transform 
non-trivially under gauge transformations Q. in the gauge group 

A^(x)^A;^(x), Q{x)^a{x)Q{x). (1.1) 

Gauge invariance tells us that we can work with any representative of the gauge orbit {A^ (x) | (x) € 
^}. A common method of choosing a unique representative is through the global maximum of a 
given gauge fixing action: 

'S'fix [A^] ^ global maximum for Q. =: Qfmr[A]. (1.2) 

The set of these representatives is free of Gribov copies and (after suitable boundary identifications) 
is called \he fundamental modular region: 

{A^{x)\a = a^^[A]]. (1.3) 

It has been known for a long time [|l|] that a physical electron state \Q) can only be obtained if the 
bare state 2(x)[0) is properly dressed by a photon cloud: 

10 = h\A]Q{x)\0), /2[a"] = /i[A]at . (1.4) 

The latter equation specifies the dressing property which ensures the gauge invariance of The 
choice of a a dressing is not unique. They may differ e.g. by a gauge invariant factor without 



spoiling the dressing property in (1.4). One particular way to construct a dressing is to use gauge 
fixing: 

h[A\ = HfmrIA] . (1.5) 

We call this definition dressing from the fundamental modular region. Under the assumption that 
there is a unique global maximum, it can be easily shown that this definition satisfies the dressing 
property: 

h[A^\ = Qfmr[A"] =^^FMR[A]a1' = h[A]a^ . 

The disadvantage of this definition is that it relies on an explicit realisation of gauge fixing. While 
the Faddeev-Popov method [§] works in perturbation theory, it fails non-perturbatively, leading to 
Green functions being in indeterminate form ^ |5|]. This is because in non-Abelian theories the 
gauge fixing condition has multiple solutions Q.{x). This was first elucidated by Gribov for the 
case of Coulomb gauge |^], and since then has been known as the Gribov problem. One might 
argue that in a purely numerical approach using lattice regularisation there is no need to resort to a 
gauge fixing condition: one might find the unique representative on the gauge orbit by seeking the 



global maximum of the gauge fixing action (1.2) using advanced algorithms. While this approach 



is conceptually correct, it is not feasible: it can be shown that locating the global maximum is 
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equivalent to solving a spin-glass problem which is beyond the scope of the numerical techniques 
currently available. 

In order to bypass the Gribov problem, we recently proposed |^] combining the alternative 
construction of the gauge invariant partition function by Parrinello and Zwanziger [||, |^ with 



the technique of gauge invariant projection initiated by Polyakov and Susskind [ ]rT| , II4 pj[ ]. This 
leads to our integral dressing construction of gauge invariant trial states, which we will describe 
below. We will also describe the "ice-limit", in which we recover the dressing from the FMR 
without having to construct the FMR explicitly. As an initial application, we employ the integral 
dressed trial states to calculate the static interquark potential. 

2. Integral dressing 

In the following, we will adopt lattice regularisation. The gauge fields (x) are represented 
by group valued links U^{x). Rather than realising the dressing property (1.4) by gauge fixing, we 
define the integral dressing |^ 

h[U] = j Q.{x) exp{K:5fix[t/"]}, (2.1) 

where K plays the role of a gauge fixing parameter. Using the gauge invariance of the Haar measure, 
we easily verify the dressing property (G G 

h[U^] = I &a a{x) exp{jc5fix[t/"^]} = I &[^G] [aG]^^ exp{jc5fix[t/"^]} = h[U]G\ 

Gauge invariant dressings for multi-quark trial states may be defined similarly, e.g. for a quark- 
antiquark trial state we make the ansatz 

Itrial) ■.= Q\y)h^^\u]{y,x) Q{x) |0), h^^\u] := j QQ. a\y)a{x) exp{jc5fix[t/'']}. (2.2) 
Let us briefly discuss the strong gauge fixing limit ^ 00. In this case, the total contribution to the 



integral in (2.2) arises solely from the domain where S^^ attains its global maximum. In this case, 



we recover the dressing from the FMR: 

oc n^^J^), \FMR) := hm |trial). (2.3) 

As an illustrative example, we consider axial dressing by choosing 

5fix[t/"] = Li^tr«73"(x) < A^,, (2.4) 

where A^^v is the number of lattice cites, and 'Nf. is the number of colours of the gauge group. The 
key observation is that the upper bound in ( |2.4| ) can be saturated exactly: letting denote the 
straight Une joining the points (xi,X2,l,X4) and (xi ,X2,X3,X4), we find HfmrW = W^g^^^i?^) in 



the strong gauge fixing limit of k" ^ 00. The integral dressing ( [2.3| ) for a quark at x and an antiquark 
aty = x + reT, then becomes 

/i(2)[[/](3;,x) = n^3W, (2.5) 
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Figure 1: Illustration of the ice-limit (left panel) and the static potential from Coulomb dress trial states 
(right panel). 



where Cxy is the straight line between x and y. Relation (2.5) and the corresponding trial state have 
an intuitive interpretation: h^^^ is just the Polyakov line between the points x and y, so the trial state 
is made gauge invariant by joining the quark and antiquark with a thin gluonic string. 

The static quark antiquark potential V{r) can be obtained by calculating the expectation value 
of the (Euclidean) time evolution operator using the FMR trial state: 



,,„;„n riarge , „n 12 ^-V Wr 



(trial|e"'"' [trial) 



|(2|trial)j^e" 



(FMR|e"^^|FMR) = W{r,T). 



(2.6) 



In the case of the axial dressing, we recover the standard approach to the static potential using 
rectangular Wilson loops W{r,t). Although we have not encountered here a Gribov problem in 
the case of axial dressing, its disadvantage is that the overlap |(2jtrial)p is poor and even vanishes 
when the lattice regulator is removed [Q, p^ ]. 



3. Dressing from the Coulomb gauge FMR - the ice-Umit 



A much better overlap is observed if Coulomb gauge fixing is employed [14|: 



1 



5fix[t/"] = iTTtrf/fW, k = \...3. (3.1) 

k,x ■'^c 

Using the definition of the integral dressing and heavy quark propagation, the expectation value 



(2.6) can be written as 



(trial le^-^^ltrial) oc J ... &x^{KS^y\U^]} 



(3.2) 



Note the order of the integration: the Q. integration is done before the link integration. With 
this ordering, we recover the Gribov problem in the large-JC limit, since for a given background 
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field U^{x), the il-path integral may be viewed as a spin-glass partition function with a metric 
U^{x). However, using the explicit gauge invariance of our approach, the order of integration can 
be exchanged, and the U integral becomes independent of Q. (for details see [^]), 

(trialje"'^^ [trial) oc (^j QQ^ J S)U^ ... exp{K-Sfix[t/]}. (3.3) 

Changing the order of integration solves the Gribov problem: the limit k" ^ oo can thus be taken 



analytically, as the dominant contributions to the integral now arise from maximising { p.l\ ) with 
respect to U rather than Q.. The solution is clearly to take U = I, i.e. 

5fix[t/] ^ max =^ Ukix,t) = 1 for f = and f = T . (3.4) 

The limit fc — > oo implies that the link fields on the time-slices t = and t = T are frozen 
to perturbative vacuum levels. We have therefore called this the ice-limit. An illustration of the 
ice-limit together with our numerical results for V (r) can be found in figure [l] - a linearly rising 
interquark potential clearly emerges. 
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